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ABSTRACT 

Context. The high star formation rates of luminous infrared galaxies (LIRGs) make them ideal places for core-collapse supernova 
(CCSN) searches. Massive star formation can often be found in coexistence with an active galactic nucleus (AGN), contributing 
jointly to the energy source of LIRGs. At radio frequencies, free from dust extinction, it is possible to detect compact components 
within the innermost LIRG nuclear regions, such as SNe and SN remnants, as well as AGN buried deep in the LIRG nuclei. 
Aims. Our study of the LIRG IC 883 aims at: (i) investigating the parsec-scale radio structure of the (circum-)nuclear regions of 
IC883; (ii) detecting at radio frequencies the two recently reported circumnuclear SNe 2010cu and 201 lhi, which were discovered 
by near-IR (NIR) adaptive optics observations of IC 883; and (iii) further investigating the nature of SN 201 lhi at NIR wavelengths. 
Methods. We used the e-EVN at 5 GHz, and e-MERLIN at 6.9 GHz, to observe contemporaneously the LIRG IC 883 at high angular- 
resolution (from tens to hundreds of milliarcsec) and with high sensitivity (< 70 /jJy), complemented by archival VLBI data at 5 and 
8.4 GHz. We also used the Gemini-North telescope to obtain late-time JHK photometry for SN 201 lhi. 

Results. The circumnuclear regions traced by e-MERLIN at 6.9 GHz have an extension of ~1 kpc, at a position angle of 130°, and 
show a striking double-sided structure, which very likely corresponds to a warped rotating ring, in agreement with previous studies. 
Our e-EVN observations at 5 GHz and complementary archival VLBI data at 5 GHz and 8.4 GHz, reveal the presence of various 
milliarcsec compact components in the nucleus of IC 883. A single compact source, an AGN candidate, dominates the emission at 
both nuclear and circumnuclear scales, as imaged with the e-EVN and e-MERLIN, respectively. The other milliarcsec components 
are very suggestive of ongoing nuclear CCSN activity. Our e-EVN observations also resulted in upper limits to the radio luminosity of 
the two SNe in IC 883 recently discovered at NIR wavelengths. We refine the classification of SN 201 lhi as a Type IIP SN according 
to our latest Gemini-North epoch from 2012, in agreement with a low-luminosity radio SN nature. We estimate a CCSN rate lower 
limit of l.ljj yr~' for the entire galaxy, based on three nuclear radio SNe and the circumnuclear SNe 2010cu and 201 lhi. 

Key words. Galaxies: starburst - Galaxies: nuclei - Galaxies: individual: IC 883 - Radio continuum: stars - Radiation mechanisms: 
non-thermal - supernovae: individual: SN2010cu, SN2011hi 



. 1 . Introduction entially concentrate (i.e ., the disk) as seen in optical and near- 

' IR (NIR) studies ( e.g., [Smith et aljri995UScoville et all I2OOOI: 

r £ ; The excess of radio continuum and infrared (IR) emission in lu- [ Modica et al.ll2012h . and a narrow one along the disk itself (Keel 

& • minous IR galaxies (LIRGs) is associated with vigorous forma- 1985) 

tion of massive stars (M ^ 8 Mo) and/or an active galactic nu- „, , . . c ,. , 

1 , . , r?; — ; — „\ ,. ', — _° Taking Lie as a measure of the rate at which massive 

cleus (AGN) (see e.g., Sanders & Mirabel 1996). Core-collapse , % m , t . . ,. ™ pv . 

L„, T r : t~! — ~ c stars are formed, results in a corresponding CCSN rate, vcn, 

supernovae (CCSNe) serve as probes of the star formation rate c , ~ _i c T „ 00 ~ , - a. • • , u 

j. . r , ■ • , •• 1 tt, /-"-< • of ~1.3yr for IC883 (assuming the empirical relation ob- 



(SFR) of massive stars, and therefore their detection in LIRGs is 
a pertinent task 



tained by Mattila & Meikle 2001). Optical spectroscopic stud- 



ies by iVeilleux et al.1 (11995) have led to the classification of 



IC 883 is a LIRG with IR luminosity (L lR - L[8-1000 fim]) the IC 883 nucl eus as a low-ionis ation narrow emission-line re 
-4.7 x 10 1 L q at a distance^ 100 Mpc (1 mas « 0.48 pc), ac- gion However, ^^rtd] dS) have recently reclassified this 
cording to | Sandersetal. | (|2003|). This LIRG is probably the re- galaxy as a composite Starbu rst-AGN, based on the revised op- 
sult of the merger of two disk galaxies (|Smifh et al.||1995|) , as tical diagnostic diagrams by iKewlev et all (120061). T he multi- 
suggested by its associated tidal tails: a diffuse one approxi- wave i e ngth study on IC883 by Modic a et al.1 (l2012h also sup- 
mately perpendicular to the region where gas and dust prefer- ports a composite nature Consequently it is expected that an 

AGN is contributing to the total Lir in IC 883, and thus the 

* e-mail: crroca@utu.fi vsn should be lower than 1.3 yr -1 . Indeed, at radio frequencies, 
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ISmith et all (Il998l) and lParra et"ai1 (l2010t) have found evidence 
of compact AGN activity through VLBI studies of IC 883. 

IC 883 is part of our program "An ALTAIR study of 
Supernovae in Luminous Infrared Galaxies", using Gemini- 
North with_jJs_Ljiser ■_Guide Star adaptive optics (AO) system 
(e.g., iKankare et al.l 120081) . During the four year survey, our 
program has yielded the discovery of two CCSN e in IC883 
within a 12-month period: SN 2010cu (24 Feb. 2010,lRvder et alj 
20101) and SN 201 l hi (11 Feb. 2011, IKankare et alJ l201lh 



Kankare et a i1 (l2012h concluded that either a Type IIP or a Iln/L 



SN provided the best fits to the NIR light curves of SNe 2010cu 
and2011hi. 

In this paper we report our radio observations towards the 
nuclear and circumnuclear regions of IC 883, as well as the most 
recent NIR epoch o btained with Gemini- North on January 31 
2012, not included in lKankare et alJd2012l) . In Section[2]we give 
details of our radio and NIR observations, Section[3]contains our 
results and discussion, and Section |4] contains the summary of 
our study. 



2. Observations 

IC 883 was observed at radio frequencies with the electronic 
Multi-Element Remotely Linked Interferometer Network (e- 
MERLIN) and with the electronic European very long base- 
line interferometry (VLBI) Network (e-EVNfl), and at NIR 
wavelengths with the Gemini-North telescope. Our observa- 
tions aimed at imaging the parsec-scale radio structure of the 
(circum-)nuclear regions of IC883, and detecting SNe 2010cu 
and 201 1 hi at radio frequencies, as well as to further investigat- 
ing the nature of SN 201 lhi at NIR wavelengths. To complement 
our radio observations, we have also analysed publicly available 
archival VLBI data. 

In Table[TJwe give some details on the observational param- 
eters of the different radio observations included in our study: e- 
MERLIN, e-EVN and additional VLBA and EVN observations. 



2.1. e-MERLIN observations and data reduction 

We carried out e-MERLIN observations of IC 883 at 6.9 GHz 
(median central frequency) on 19 March 2011. These were 
Director's Discretionary Time (DDT) observations within the 
commissioning phase of e-MERLIN, which included the fol- 
lowing 25 m diameter antennas: Mark II, Defford, Knockin, 
Damhall and Pickmere. The observations lasted ~ 24 hr, from 
which approximately 20 were spent on target and 10 hr were ulti- 
mately usable after editing. Four sub-bands (512 channels each) 
with dual polarisation were used, accounting for a total band- 
width of 5 1 2 MHz. 

We analysed the data within the NRAO Astronomical Image 
Processing System (aips). 3C286 set the absolute flux density 
scale following an iterative process. The flux estimated for 
3C286 in the shortest baseline (Mark II-Pickmere), at the centre 
of the different sub-bands, was on average 5.56 Jy. We then cali- 
brated the amplitude of DA 193, which is a bright and unresolved 
source as seen by all e-MERLIN baselines, resulting in an aver- 
age flux density of 3.72 Jy. Finally, we used DA 193 to set the 
flux density of the phase reference source, J 1324+3622, which 
resulted in an average flux density of 68.25 mJy. We performed 



1 The development of e-VLBI within the EVN has been made 
possible via EXPReS project funded by the EC FP6 1ST Integrated 
Infrastructure Initiative contract # 026642 - with a goal to achieve 1 
Gbps e-VLBI real-time data transfer and correlation. 



a series of phase-only self-calibration iterations of the phase ref- 
erence source, before phase-calibrating the target source. To ac- 
count for the correlation offset from the position of the strongest 
source in the field, we used the task UVSUB in aips. To improve 
the sensitivity, we performed phase-only self-calibration itera- 
tions on the target source. We achieved a thermal rms noise in the 
IC 883 map of 44 /iJy/beam, for a beam size of 165x88 mas at 
- 1 1 .6° (see left panel in Figure[TJ, obtained with natural weight- 
ing and ROBUST=l within aips. 

2.2. e-EVN observations and data reduction 

We observed IC883 on 23 March 2011 (ToO project: RR006, 
PL: Romero-Canizales) at 5 GHz, with the e-EVN, which in- 
cluded the following antennas (diameter, location): Effelsberg 
(100 m, Germany), Mark II (25 m, UK), Medicina (32 m, Italy), 
Onsala (25 m, Sweden), Torun (32 m, Poland), Westerbork array 
(14x25 m, NL) and Yebes (40 m, Spain). 

RR006 was a 2 hr experiment (~ 1 .3 hr, total time on source), 
recorded at 1024 Mbps using eight sub-bands, each of 16 MHz 
and dual polarisation. The data were correlated at the EVN 
MklV Data Processor at JIVE with an averaging time of 2 s. The 
point-like source Jl 159+2914 (2.45 Jy at 5 GHz) was used as 
a fringe finder and bandpass calibrator. J1317+3425 (0.35 Jy at 
5 GHz), at ~0.7° angular distance from IC 883, served as phase 
reference source. 3.4 min scans on IC883 were alternated with 
1.2 min scans on J1317+3425. The data were ana lysed with aips. 
We a lso used the Caltech program DIFMAP (Shepherd et al. 
1995) to image the calibrators and to assess the performance of 
each antenna. 

In an attempt to detect a radio counterpart for SN 201 lhi, its 
coordinates derived from the NIR images were used for pointing 
and correlation. A local rms of ~30 yuJy/beam was achieved at 
the pointing centre. Since no radio source was detected above 
3cr at that position, we corrected the visibility data with the task 
UVSUB in aips, as we did for the e-MERLIN data, to account 
for the correlation offset from the position of the strongest source 
in the field, i.e., to the radio nucleus of IC883. The shift of 
~0.8 arcsec is well within the primary beam of each antenna in 
the used array. 

For the imaging process within aips we used different 
weighting schemes in order to test the reliability of our results. 
Due to the scarce uv-coverage, the use of natural weighting tends 
to give rise to secondary side lobes at a 60 per cent level of 
the main lobe, which can produce artifacts in the reconstruction 
imaging process that can be confused with putative sources. Our 
final image was made using uniform weighting and ROBUST= 1 . 
We achieved a final thermal rms in the map of 66.2 /vJy/beam, 
for a beam size of 9.20x6.36 mas at -76.1° (see Figure Q] right 
panel). 

2.3. Archival VLBI data 

We have analysed Very Long Baseline Array (VLBA) data from 
projects BN026 and BN027 at 5 GHz, and BC196 at 8.4 GHz. 
The VLBA observations under pr ogram BC196 were processed 
with the DiFX VLBI correlator (iDeller et al.l 1201 ll) . We also 
report EVN observations at 5 GHz made on June 15th 2006 
(project: EP055, PL: R.Parra); see Table[TJ 

In project EP055, IC 883 was observed in four scans spread 
over 24 hr to optimize the uv-coverage for a total observation 
time on source of ~1.3hr. The array used comprised 10 EVN 
stations: Lovell (76 m, UK), Cambridge (32 m, UK), Westerbork 
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Table 1. Observational parameters of the radio observations. 



Epoch 


Project 


Array 


Observing 


V 


Phase 


Py 


Weighting 


Convolving beam 


mis 


label 






date 


(GHz) 


calibrator 


( Jy/beam) 




(mas) 


( //Jy/beam) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


EM 




e-MERLIN 


2011-03-19 


6.9 


J1324+3622 


0.07 


N,l 


165.23x88.35 at 11.6° 


44.0 


VI 


BN026 


VLBA 


2004-08-13 


5.0 


J1317+3425 


0.25 


N,0 


3.71x1.47 at -0.1° 


84.7 


V2 


BN027 


VLBA 


2004-09-20 


5.0 


J1310+3220 


1.41 


N,l 


3.48x1.65 at 17.6° 


119.0 


V3 


BN027 


VLBA 


2005-07-11 


5.0 


J13 10+3220 


1.19 


N,2 


3.43x1.78 at 11.5° 


107.0 


V4 


EP055 


EVN 


2006-06-15 


5.0 


J2333+3901 


0.81 


U,l 


7.99x6.00 at -81.3° 


162.0 


V5 


RR006 


e-EVN 


2011-03-23 


5.0 


J1317+3425 


0.35 


U,l 


9.20x6.36 at -76.1° 


66.2 


V6 


BC196 


VLBA 


2011-05-15 


8.4 


J2330+3348 


0.33 


N,0 


2.41x1.17 at -8.3° 


174.0 



Notes. Phase reference sources used in each observation are listed along with their associated peak intensities (columns 6 and 7). For the different 
IC 883 maps, column 8 lists the weighting parameters used in the imaging process: N=natural, U=uniform, and the number corresponds to the 
value of the ROBUST parameter). Column 9 lists the resolution, and column 10 the attained rms noise. 



array (14x25 m, NL), Effelsberg (100m, Germany), Medicina 
(32 m, Italy), Noto (32 m, Italy), Onsala (25 m, Sweden), Torun 
(32 m, Poland), Urumqi (25 m, China) and Shanghai (25 m, 
China). The observations were performed in a 1024 Mbps dual 
polarization mode. The data were correlated at JIVE in the 
Netherlands with an averaging time of 2 s. 

All the additional VLBI data were edited and reduced using 
standard procedures within aips. The phase reference calibrator 
used in project BN027 is slightly resolved at VLBA scales, and 
we have thus corrected the fringe solutions from its structure. 
We did not apply any self-calibration to obtain the final images 
of IC883. In the case of project EP055, the data from the two 
Chinese antennas (Ur and Sh) as well as from Cm could not be 
properly calibrated, so we omitted these antennas in the imaging 
process. 

As in the case of our e-EVN epoch, we have tested the re- 
liability of the sources in the different VLBI maps, by using 
different weighting schemes. We used the atps task UVSUB to 
apply shifts in a and 5 to the visibility data on all the addi- 
tional archival VLBI epochs. This allowed us to have all the 
maps centred at the same position (a(J2000) = 13 /l 20'"35 s 3184, 
(5(J2000) = 34°08'22'.'352). Since the resulting images have a 
different convolving beam and were thus mapped with different 
pixel sizes, we changed their geometry to match that of our e- 
EVN epoch (pixel size of 0.5 mas), by means of the task OHGEO 
within aips. 

2.4. Gemini-North observations 

Our Gemini-North observations of IC 883 were conducted with 
the ALTAIR Laser Guide Star AO system on the Near-InfraRed 
Imager (NIRI). The LIRG was observed in 1 1 separate epochs 
between 2008 April 15.5 UT and 2012 January 31.6 UT, cov- 
ering a period of 3.8 yr. Here we report the results from the 
most recent epoch of JHK imaging obtained in 2012 under pro- 
gram GN- 2012A-Q-56 (PL: S.R yder); the previous epochs are 
reported in Kankar eet al.1 (120121) . For details on the data reduc- 
tion, calibration, image subtraction process and photometry of 
the ALTAIR/NIRI da ta on SN 2010cu and SN 201 lhi we refer to 
iKankareet all (120121) . 

Image subtra ction based on the ISIS 2.2 dAlard & Luptonl 
ll998HAlardll2000l) package was used to compare the 2012 JHK 
images of IC883 to the reference data obtained on 2010 May 
4-5 UT. This revealed SN 201 lhi to have faded below the detec- 
tion limit of the Gemini imag es. Similar to deriv ing the /-band 
upper limit for SN 2 01 Ocu in IK ankare et we used the 

QUBA package (see IValenti et al.ll201 II) running standard iraf 



tasks, to measure the faintest detectable sources in the 2010 ref- 
erence image with poorer quality AO correction compared to the 
latest 2012 epoch. Taking into account the noise increasing by a 
factor of V2 due to the image subtraction process, conservative 
5<x upper limits of roughly mj > 19.7 mag, mn > 19.9 mag, and 
m.K > 19.9 mag were derived. 

3. Results and discussion 

The images of the nuclear and circumnuclear regions of IC 883 
as observed with the e-EVN and e-MERLIN are shown in Figure 
[TJ In Figure [2] we show the images resulting from a total of six 
VLBI epochs, covering a period of seven years. In Table [2] we 
report the estimated parameters for the > 5<x sources detected in 
all the images. 

Our new Gemini-North epoch has yielded new JHK magni- 
tude limits for SN 201 lhi. In Figure we show the SN 201 lhi 
template light curve fits from iKankare et al.l (1201 2l) . including 
the new limits obtained from our 2012 observations for compar- 
ison. 

We note that the discovery of two SNe (2010cu and 201 lhi) 
in a 3.8 yr period results in a vsn ~ 0.5^^3 yr~' (using the up- 
per and lower Poisson l<x uncertainties for two events given by 
lGehrelsl fl986). which agrees within the uncertainties with the 
IR luminosity-based CCSN rate of the galaxy. However, this es- 
timate does not include any correction factors nor control time 
considerations, and should thus be regarded as a lower limit. 

3.1. The nature of SNe 201 Ocu and 201 1hi 

In IKankare et al] ( |2012|) we investigated the nature of SNe 
201 Ocu and 201 lhi by means of five epochs of NIR photomet- 
ric observations prior to 2012. Based on the comparison of the 
light curves of SN 2010cu and SN 201 lhi to the CCSN tem- 
plates presented in Mat tila & Meikld (j2001) and the well sam- 
pled NIR light curves of the canonical Type IIP SN 1999em 
dKrisciunas et alj 12009), it was concluded that t he best fit for 
both S Ne was either a Type IIP or a Iln/L SN (see Kankar e~et al.l 
120121) . The ordinary, linearly declining template fit was previ- 
ously excluded because such fit suggested an extremely early 
discovery and an unrealistically long rise time to the light curve 
peak. A NIR-brigh t Type Hn/L SN w ith ~0mag of extinction 
was favoured in Kankar e et al.l (1201 2|) as the type of SN 201 lhi 
due to the low (< 2 mag) total line-of-sight extinction implied 
by the J-K colour map of the galaxy. However, a Type IIP light 
curve with higher ex tinction could also pro vide a good fit to the 
data, and as noted in Kankar e~et al.l d2012l) . higher localized ex- 
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Fig. 1. IC 883 contour images at a median central frequency of 6.9 GHz obtained with e-MERLIN (left), and at 5 GHz with the e- 
EVN (right), in observations carried out in March 201 1. The rms noise in the two images is 44 and 66 yuJy/beam, respectively, and 
the convolving beam in each case is indicated the upper left corner of the maps. Dashed contours represent -3cr levels. Note that the 
brightest component at e-MERLIN scales is dominated by a compact source seen at mas scales with the e-EVN. Th e crosses in the 
left p anel indicate the coordinates and positional errors of SNe 2010cu and 201 lhi, according to the values reported in lKankare et al.l 
(2012). 



tinctions are possible, as the colour map is tracing only the large 
scale structures of the LIRG. 

In Figure |3]we present the template y 2 fits o f the SN2011hi 
light curves adopting the ICalzetti et al.l (|2000) reddening law 
and including the new upper limits. However, these limits are 
not i ncluded in the fit, i. e., the values reported in the Table 
2 of iKankare et alj d2012l) are unchanged. Based on the new 
data, the slowly declining (NIR-bright Type Iln/L) SN tem- 
plate is now excluded due to the fading of the SN. In addi- 
tion to the SN 1999em template, we have now also compared 
the NIR fight curve of SN 201 lhi w ith that of other Type IIP 
SNe: 2003hn dKrisciunas et all2009h and 2004et dMaguire et alJ 
12011 . The results of the best fits are similar to those obtained 
with SN 1999em, but having larger^ 2 values. The V-band ex- 
tinctions from the three different Type IIP SNe used for compar- 
ison are consistent within 0.5 mag for the adopted extinction law. 
We conclude that SN 201 lhi was in fact a slightly over-luminous 
(1-2 mag brighter than SN 1999em) Type IIP SN discovered 
roughly a month after the explosion, with 5 or 7 mag of host 
alaxy extinc tion in the V-band ad opting either the lCardelli et al.l 



galaxy 
(1989) 



or the Calzetti et al. (2000) reddening law, respectively. 

There are no radio sources detected with the e-EVN above 
3cr at/around the reported positions for either SN2010cu, nor 
SN 201 lhi. Therefore, the 3<x upper limits at 5 GHz correspond 
to 2.4 x 10 27 erg s" 1 Hz 1 for both SNe. The SNe are also not 
detected in the archival 8.4 GHz VLBA observations (project 
BC196) carried out two months after our e-EVN observations. 
Source B2a detected with e-MERLIN at 6.9 GHz is positionally 
consistent within the uncertainties with SN2010cu. However, 
B2a is not a compact source, and we thus discard the possibility 
of it being a SN. 



SN 2010cu, has an estimated age of 14-17 months at the time 
of the radio observations made in 2011. Its non-detection at 5, 
6.9 and 8.4 GHz is an indication that: (i) the radio emission at 
these frequencies was well into the optically thin decline phase 
of its evolution (e.g., a Type Ib/c SN), now having a luminos- 
ity well below 2.4 x 10 27 ergs" 1 Hz" 1 ; and/or (ii) SN2010cu is 
an intrinsically low-luminosity Type II radio supernova (e.g., a 
Type IIP such as SN 1999em) and thus more sensitive observa- 
tions would be required for its detection. 

SN 201 lhi was also not detected in our radio observations 
made 2-5 months after its explosion, and this is consistent 
with the fact that type IIP SNe are low-luminosity radio emit- 
ters. In fact, other we ll known type IIP SNe like SNe 1999em , 
2004di. and 2004et (jPoolev et alJ l200l iBeswick et ail l2005b 
Marti- Vid al et alJl2007l respectively) have peak radio luminosi- 
ties b elow 10 26 erg s" 1 Hz" 1 (see e.g., figure 2 in lChevalier et al.l 

2006) . At NIR wavelengths SN 201 lhi was only slightly brighter 
than SN 1999em; if its radio emission was likewise only slightly 
stronger, we can expect that its radio luminosity was well below 
the 3cr detection threshold (2.4 x 10 27 erg s" 1 Hz" 1 ) of our most 
sensitive radio observations carried out with the e-EVN. We note 
that the non-detection of SNe 2010cu and 201 lhi at radio wave- 
lengths does not imply that no radio emission has been produced 
by these NIR detected SNe. 

To gain a better understanding of SNe in LIRGs, it is desir- 
able to have a radio monitoring of the NIR-detected SNe includ- 
ing observations at different frequencies over many epochs. This 
strategy would help to detect either a fast or a slowly-evolving 
SN, w hich can also help to constrain their type (Stockdale et al.l 

2007) . However, not all CCSNe will be bright radio emitters. 
This depends on the interaction between the SN ejecta and the 
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Fig. 2. IC 883 contour maps of the VLBI observations listed in Table Q] All the maps are centred on a(J2000) = 13''20'"35 s :3184, 
(5(J2000) = 34°08'22'.'352. The crosses indicate the positions from the detection epoch of the different 5<x sources reported in Table 
[2) Al and A2 from experiment V5; A3 from experiment V6; and A4 from experiment V4. 
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circumstellar medium, which give s rise to the sync hrotron emis- 
sion at radio frequencies (see e.g.. IChevalierll 1 9 82b . 



3.2. The circumnuclear regions of IC 883: a jet or a disk? 

The quest to detect radio emission from SNe 2010cu and 201 1 hi 
has brought the opportunity to unveil different facets of IC 883. 



In the following we discuss the circumnuclear structure of 
IC 883 and the nature of its nuclear regions. 

Our e-MERLIN observations reveal a double-sided structure 
with an approximate projected extent of 2 arcsec (~ 1 kpc) by 
0.3 arcsec (« 150 pc), at a position angle of 130° (Figure[TJ. This 
inclinatio n matches that prev i ously found in o ther radio stud- 
ies (e.g.. ICondon et alJ 1 199 it iParra et aO 120 10). However, the 
double-sided structure in IC 883 is revealed for the first time at 
6.9 GHz with our e-MERLIN observations (see also the 1 .4 GHz 
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Table 2. Estimated parameters of the detected radio emitting sources in the (circum-)nuclear regions of IC 883. 



Source 


Epoch 


Aa'(J2000) 


A<5"(J2000) 


Py 


Sy 


L v 




©M x & m 


logr B 


name 


label 






( mjy/beam) 


(mJy) 


(10 27 ergs- 


'Hz" 1 ) 


(mas 2 ) 


(K) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


A 


EM 


35.319 (<1) 


22.35 (1) 


4.89 ± 0.25 


6.07 ±0.31 


72.65 ± 


3.67 


104.8 x 30.4 


4.15 ± 0.02 


Bla 


EM 


35.332 (1) 


22.15 (3) 


1.04 ± 0.07 


2.07 ±0.11 


24.75 ± 


1.34 


138.1 x 97.1 


3.45 ± 0.02 


Bib 


EM 


35.304 (1) 


22.63 (6) 


0.59 ± 0.05 


1.75 ±0.10 


20.97 ± 


1.17 


191.3 x 141.6 


3.09 ± 0.02 


B2a 


EM 


35.356 (1) 


21.74 (3) 


1.03 ± 0.07 


2.54 ± 0.13 


30.34 ± 


1.61 


224.0 x 90.5 


3.11 ±0.02 


B2b 


EM 


35.281 (2) 


23.05 (11) 


0.31 ±0.05 


0.97 ± 0.07 


11.55 ± 


0.78 


295.9 x 109.8 


2.45 ± 0.03 


Al 


VI 


31.3180 (-0.0) 


22.352 (0.1) 


1.68 ±0.12 


1.94 ±0.13 


23.19 ± 


1.54 


1.2 x 0.6 


7.81 ±0.03 




V2 


31.3179 (0.1) 


22.352 (-0.0) 


1.53 ±0.14 


1.53 ±0.14 


18.33 ± 


1.69 


0.6 x ■ ■ ■ 


6.79 ± 0.04 




V3 


35.3181 (0.1) 


22.349(0.1) 


1.95 ±0.15 


2.65 ± 0.17 


31.73 ± 


2.04 


1.8 x 0.3 


7.59 ± 0.03 




V4 


35.3180 (0.6) 


22.351 (-0.0) 


1.12 ±0.17 


1.21 ±0.17 


14.49 ± 


2.07 


13.0 x ■■■ 


5.97 ± 0.06 




V5 


35.3179 (0.1) 


22.352 (-0.0) 


3.82 ± 0.20 


3.82 ± 0.20 


45.67 ± 


2.42 


0.8 x ■■■ 


6.35 ± 0.02 




V6 


35.3178 (-0.0) 


22.351 (0.1) 


2.85 ± 0.22 


4.38 ± 0.28 


52.34 ± 


3.34 


1.5 x 0.4 


7.55 ± 0.03 


A2 


V5 


35.3190 (0.7) 


22.317 (0.2) 


0.42 ± 0.07 


0.56 ± 0.07 


6.66 ± 


0.86 


7.0 x 2.0 


5.75 ± 0.06 


A3 


V6 


35.3176 (-0.0) 


22.354 (0.2) 


0.97 ±0.18 


1.55 ±0.19 


18.54 ± 


2.28 


1.4 x 1.1 


7.15 ±0.05 


A4 


V4 


35.3181 (0.7) 


22.371 (0.1) 


0.89 ±0.17 


1.27 ± 0.17 


15.25 ± 


2.08 


8.1 x 0.6 


5.98 ± 0.06 




V5 


35.3187 (0.9) 


22.366 (0.2) 


0.33 ± 0.07 


0.67 ± 0.07 


8.06 ± 


0.89 


14.0 x 1.4 


5.23 ± 0.05 



Notes. Columns: (1) Source names corresponding to those in Figures [TJand [2] (2) Epoch label corresponding to those specified in Table[TJ Note 
that epoch EM and V5 are quasi-simultaneous. (3^1) Coordinates given with respect to o-(72000) = 13' ! 20'"0O!0 and (5(72000) = 34°08'0O70. 
Position uncertainties in mas, within parentheses, are given by FWHM/(2 x SNR), where SNR is the signal to noise ratio, and FWHM was taken 
as the projection of the beam major axis on both a and S axes. (5) Peak intensity. We estimated the uncertainties by adding in quadrature the rms 
noise in the map (column 10 in Table[T} plus a 5% uncertainty in the point source flux density calibration. (6) Flux density. (7) Monochromatic 
luminosity at the observed frequency (column 5 in Table[T). (8) Deconvolved major and minor axes, obtained by fitting a Gaussian to each source. 
(9) Brightness temperature. When no deconvolved sizes could be obtained, we used instead the solid angle subtended by the synthesised beam for 
calculating lower limits of J B , i- e -> &s = 7r(41og2)~' (FWHM M x FWHM m ), with FWHM M and FWHM m , the major and minor synthesised beam 
fitted FWHM (column 9 in Table[T}, respectively. 



MERLIN map fro m IClemens & Alexanderil2004l) VLA obser- 
vations (e.g., from|Condon et al Jl 1 99 ll iParra et al.ll20 1 Ol) do not 
have enough spatial resolution to image the circumnuclear re- 
gions with similar level of detail. By convolving our e-MERLIN 
image with a beam of size (X'4, which corresponds to VLA A- 
configuration observations at 5 GHz, we were able to recover 
a struc ture similar to that shown in figure 2 from IParra et aTl 
(I2010h . 



At first glance, the structure seen with e-MERLIN resembles 
that of a two-sided jet of an AGN, whose restarting activity ex- 
plains the presence of two sets of condensations (Bla-Blb and 
B2a-B2b; see left panel in Figure [TJ, at either side of the core 
A. The larger flux density of the approaching components in the 
jet, Bla and B2a in this case, with respect to the flux density of 
the receding components, Bib and B2b, could be explained by 
Doppler boosting effects. 

There are however strong arguments against a jet interpre- 
tation: (i) the approaching components are at shorter distances 
from the core than the receding components, thus implying that 
their apparent velocities are smaller; (ii) the putative jet struc- 
ture is placed on top of the disk of the galaxy, and this is difficult 
to explain, unless we consider that the directions of jet and disk 
coincide due to a projection effect. 

Diffe rent studies of the molecul a r and atomic gas towards 
IC883 dDownes & Solomon! Il998t iBrvant & Scoviiiel 1 19991: 
IClemens & Alexander 2004), give evidence of the rotation of the 
disk from South-East towards North- West. This situation would 
imply that the putative jet is perpendicular to the rotation axis, 
and hence no projecti on effect could expla in the superposition 
of the jet on the disk. Schm itt et alj d2002) studied the orienta- 
tion of jets with respect to the rotation axis in a sample of radio 
galaxies and found that large misalignments (< 77°) may occur 
due to warping mechanisms in the black hole's accretion disk. 
In the case of IC 883 the misalignment would be close to 90°, so 



a jet scenario for the structure revealed with e-MERLIN is even 
less viable. 

If we consider instead a rotating disk (or ring) scenario, com- 
ponent A would be at the centre, while Bla, B2a (approach- 
ing components) and Bib, B2b (receding components) would be 
continuum radio components of the disk itself. We note that the 
line joining component Bla with Bib through A, is not the same 
as the line joining B2a with B2b, which gives evidence of the 
disk/ring being warped (in ag reement with previous studies, e.g., 
IClemens & Alexander! 12004). probably due to the interaction 
with the ambient material. We note that the structure unveiled 
by e-MERLIN has a width of approximately 150pc. Whilst such 
width might appear too narrow for representing a disk, thin disks 
of a few hundreds of parsec s are not uncommo n for edge-on 
galaxies at radio frequencies (Du mke et all 1995). Furthermore, 
e-MERLIN resolves out some of the extended emission which is 
otherwise visible at lower resolution, for instance with the VLA 
(ICondon et al.ll991l:|Parra et al.l2010l) . and hence, the size of the 
structure mapped with e-MERLIN is compatible with a disk/ring 
structure. 

Thus, whilst we cannot completely rule out the possibility 
that the structure mapped by e-MERLIN is a jet, there is strong 
evidence against it, and evidence favouring a rotating disk/ring 
nature instead. 

Further e-MERLIN observations at different frequencies 
could be used to probe the radio lifetime of the condensations at 
either side of the core A. Observations optimized for polarime- 
try studies would yield information on the configuration of the 
large-scale magnetic field, which could probe the existence of 
shocked material as expected in a jet scenario, hence helping to 
arrive to a firm conclusion. 
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3.3. The ongoing nuclear activity in IC 883 

Our e-EVN observations (see right panel in Figure [1}, made si- 
multaneous to our e-MERLIN observations, reveal the presence 
of at least two compact sources (labelled as Al and A2) above 
5<r in the innermost nuclear regions of IC 883. The high bright- 
ness temperature of sources Al and A2 indicates a non-thermal 
origin for both of them. 

The brightest component in our e-MERLIN image (A) has a 
peak intensity similar to that of A 1. By adding up the flux densi- 
ties of components Al and A2 in our e-EVN map, we recover a 
flux density of 4.38 ± 0.27 mJy, which is still lower than both the 
peak intensity and the flux density of component A (see Table 
|2J. We find that Al has a strong contribution to the flux density 
at larger scales (although seen at slightly different frequencies), 
and furthermore, has a position highly coincident with compo- 
nent A, despite the use of different phase reference sources and 
different spatial resolutions. Hence, Al is very likely the main 
powering source for the radio emission of component A. 

From our e-EVN observations alone we cannot characterise 
the nature of sources Al and A2, nor investigate a possible rela- 
tion between components Al and A. We have therefore searched 
for additional, publicly available, archival VLBI data (see Table 
0} that could be of use to investigate these issues. 

Component Al from our e-EVN observations is detected in 
all the additional VLBI observations (see Figure |2]l and appears 
to fluctuate in flux density at 5 GHz, varying for at least a factor 
of three between the minimum and maximum value displayed 
within the seven years comprised by the different VLBI epochs 
(see Table|2|. Component A2 is not detected in any other epoch. 
At 8.4 GHz (experiment BC196) we detect a 5<x source (A3) 
North- West of component Al. The resolution of our e-EVN ob- 
servations (RR006) does not allow us to detect component A3, 
but the source is also not detected in any of the other VLBA 
epochs at 5 GHz, which have a resolution similar to that of ex- 
periment BC196. A 3cr component placed North-East from com- 
ponent Al, is detected as a 5 <x source in experiment EP055 and 
is labelled as A4. This source has no counterpart at 8.4 GHz in 
the observations under project BC196, indicating that at that fre- 
quency the emission from A4 has become optically thin, as ex- 
pected in a SN scenario. We thus have a total of three transient 
sources (A2, A3 and A4) in the nucleus of IC 883, in addition 
to the more luminous, variable, component Al, which has re- 
mained visible for at least seven years. 

3.4. The nature of the compact sources in the IC 883 nucleus 

The milliarcsec sources detected in the innermost regions of 
IC 883 give evidence of its ongoing nuclear activity. A variable 
source, Al, has been detected in six different VLBI epochs (see 
Tableland Figure |2]i throughout seven years, and appears to be 
the main powering source of the emission at larger scales im- 
aged with e-MERLIN. Three transient sources have also been 
detected. All the sources have sizes and brightness temperatures 
indicating a non-thermal origin. The transient sources are most 
probably SNe and/or SN remnants (SNR), whereas for compo- 
nent Al the nature is less clear. To investigate the nature of 
sources A1-A4, we compare their luminosities with those re- 
ported in the literature for the non-thermal sources (SNe and 
SNRs) at 5 GHz in the starburst galaxy M82, the LIRG Arp299- 
A and in the ULIRG Arp 220 (Figure |4} . 

From Figure|4]we see that the nuclear components of IC 883, 
although being less numerous, are in general more luminous than 
the components identified as SNe and SNRs in the nuclear re- 



gions of other galaxies with high star formation rates. The life- 
time of the transient sources in IC 883 ranges from less than a 
few months up to a few years. It is thus very likely that these 
sources are a mixture of SNe and SNRs, therefore indicating a 
recent starburst activity in the nucleus of IC 883. We lack spec- 
tral information that could have helped to differentiate between 
a SN or a SNR origin for each of these transient sources. 

Component Al is the most luminous of all the sources (see 
histogram in Figure|4]i. A SN or a SNR origin are not very likely, 
based on its flux density fluctuations (see columns 5 and 6 in 
Tabled, its longevity, and its relatively flat two-point spectral in- 
dex between 5 and 8.4GHz, a = 0.3 (S v ~ v"), calculated from 
the flux densities obtained in experiments RR006 and BC196. 
Therefore Al is very likely an AGN. 

Further evidence for an obscured AGN in IC 883 can be in- 
ferred from the rati o of radio to hard X-ray luminosity (Rx = 
vL v (5 GHz)/L HX ; iTerashima & Wilson! 120031) . Ilwasawa et all 
(1201 ll) report a hard X-ray luminosity (Lhx = L[2-10keV]) 
of 6.4 x 10 40 ergs -1 for IC883. Considering the luminosity at 
5 GHz of Al observed with the e-EVN (see Table we calcu- 
late Sx ~ 3.6 x 10~ 3 , which is consistent with a low-luminosity 
AGN (LLAGN) or a normal AGN in a radio galaxy, as inferred 
from figure 4 of ITerashima & Wilsonl (120031) . 

By comparing the maximum variability (a factor of ~3) dis- 
played by Al in our observations, with the largest inter- year 
variability displayed by typical LLAGNs dNagar et al.l2002l) . we 
find that Al would be placed at the high end of the LLAGN be- 
haviour. Thus, we can only conclude that Al is a candidate for 
either a normal AGN or a low-luminosity one. 

Component A4 whose elongation could indicate a sort of 
outflow, has however no counterpart at 8.4 GHz (experiment 
BC196). Future VLBI observations at different frequencies 
should allow an unequivocal characterisation of the AGN can- 
didate in IC 883, as well as possibly discovering further nuclear 
components. 

If we consider sources A2, A3 and A4 as having their ori- 
gin in recent star formation activity, we would have three CCSN 
events that have occurred in the period between June 2006 and 
May 201 1, i.e., five years. We can thus estimate a CCSN rate for 
the IC883 nuclear region of w 0.6^ 3 yr , using the u pper and 
lower Poisson lcr uncertainties for three events given bv lGehrelsl 
(119861) . Combining the nuclear CCSN rate with the CCSN rate 
inferred from the detection of SNe 2010cu and 2011 hi within 
3.8 yr in the circumnuclear regions of IC 883, we obtain a CCSN 
rate for the entire galaxy of w 1.1*™ yr" . These rates are con- 
sidered as lower limits since our irregular temporal sampling 
both at radio and at NIR wavelengths might have missed rapidly 
evolving SNe, and also because our 5<x radio luminosity detec- 
tion threshold (4 x 10 27 erg s _1 Hz -1 ) prevents us from detecting 
Type IIP and lib SNe, and even some Type ILL. We note that the 
total rate estimate is in agreement with the IR luminosity-based 
CCSN rate estimate (see Section[TJ. It is thus of great interest to 
estimate the contribution of a putative AGN to the IR luminosity 
in order to obtain a more reliable estimate of the expected CCSN 
rate. 

3.5. The infrared SED of IC883: estimating the contribution 
of an AGN 

As mentioned in Section Q] IC 883 has been classified as a 
composite s tarburst-AGN, based on optical diagnostic diagrams 
dYuan et al.ll2010t) . 

We have modelled the spectral energy distribution (SED) 
of IC 883 (see Figure |5]l to estimate the contribution of an 
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Fig. 4. Luminosity distribution at 5 GHz of SNe and SNRs in 
M82, Arp299-A, Arp220. All the IC883 nuclear sources have 
been added for comparison purposes. Each bin has a width of 
logio(L5GHz) = 0.2 with LsGHz in units of ergs -1 Hz -1 . The SN 



and S NR l uminosities have bee n retrieved from Fenech et al. 
(2008f) and iMarchili et all (l2010h for M82. iBondi et al.1 d2012h 
for Arp299-A and from |Batejateta!J (|2011j) for Arp220. The 
most luminous M82 source corresponds to SN2008iz. In the 
case of Arp220 we have included SNe and SNRs from both 
West and East nuclei. The luminosities of the IC 883 compo- 
nents were taken from Table [2] from experiment RR006 for Al 
and A2, from experiment BC196 for A3, and from experiment 
EP055 for A4. We have assumed a spectral index of -0.7 to ob- 
tain the luminosity at 5 GHz of component A3. Note that Al is 
the brightest source, and with our most sensitive VLBI image 
(RR006) we are only sensitive to sources above a 5<x luminosity 
limit of 4 x 10 27 erg s _1 Hz -1 , indicated by a dashed line in the 
histogram. 



AGN to the total IR luminosity and to provide an estimate 
for the CCSN rate in this galaxy. For this we used a grid 
of AGN dusty torus models that have been c omp uted with 
the method of lEfstathiou & Rowan-Robinsonl dl995h and a 
grid of st arburst models that h ave been computed with the 
method of Efstathiou et al. (2000), but with a revised dust model 
(Efstathiou & Siebenmorgen 2009). The spatial resolution of the 
mid- and far-IR data used for the modelling is not high enough 
to spatially separate the different components. Therefore to dis- 
tinguish the contribution of the AGN torus and the starburst we 
have fitted a model combining both an AGN and a starburst. 

For the AGN, we considered a grid of tapered disc mod- 
els co mputed with the method of Efstat hiou & Rowan-Robinsonl 
(fl99l . In this grid of models we consider four discrete values 
for the equatorial 1000 A optical depth (500, 750, 1000, 1250), 
three values for the ratio of outer to inner disc radii (20, 60, 100) 
and three values for the opening angle of the disc (30°, 45° and 
60°). The best fit parameters which resulted in the minimum^- 2 
when we combine the starburst and the AGN torus models, are 
an opening angle of 60° for the torus, an equatorial optical depth 
of 1250 at 1000 A, and a ratio of outer to inner radius of 60. The 
AGN torus is viewed close to edge-on due to the IC 883 geome- 
try; thus, its apparent contribution to the total infrared emission 




Log (A//izm) 

Fig. 5. The infrared SED of IC 883 inc luding 2MASS and IRAS 
data, as well as the SCUBA point from lDunne et al.1 (l2000h . The 
red-dashed line is our starburst model fit, and the blue-dotted 
line is the model fit for an AGN upper limit contribution. The 
black-solid line is the model fit including both components. 



is insignificant (probably <10 per cent of the starburst luminos- 
ity). Therefore the starburst clearly dominates the observed lu- 
minosity of the entire galaxy. We note that due to viewing the 
torus almost edge-on and due to its extremely high optical depth, 
the intrinsic luminosity of the AGN itself could be a factor of 
~30 times higher than the observed one. This enhan cing factor 
is rela ted to the anisotropy of the torus emission (see Efst athioul 
2006, for a definition), and implies that the intrinsic luminosity 
of the AGN could be about three times higher than the observed 
starburst luminosity. 

Our model incorporates the s tellar population synthesis 
model of Bruzual & Charlotl (120031) . which makes a prediction 
of the SN rate, vsn(0> at a time t after star formation in an in- 
stantaneous burst. The starburst model of Efstathiou et al. (2000) 
predicts the spectrum of this instantaneous burst at time t and 
assumes a star formation history for the starburst. We assume 
that the SFR in IC 883 declines exponentially with an e-folding 
time of 20 Myr, and that the age of the starburst is 55 Myr. The 
initial optical depth of the molecular clouds that constitute the 
starburst is assumed to be 100. The SFR averaged over the du- 
ration of the starburst is 185M yr _1 . It is also possible to cal- 
culate self-consistently SN rate at different stages in the evolu- 
tion of a starburst by convolving the star formation history with 
vsn(0- This yields a vsn of 1.1 yr~' which is relevant for the en- 
tire galaxy and can therefore be compared with the observed val- 
ues. The combination of the CCSN rate for the nuclear regions 
based on three radio SNe, and the CCSN rate for the circumnu- 
clear regions based on the discovery of SNe 2010cu and 201 lhi, 
resulted in a lower limit for the CCSN rate of a; 1.1 yr~ 1 for 
the entire galaxy (see Section l3~4l >. This CCSN rate estimate is 
in good agreement with the predicted CCSN rate for the entire 
IC 883 galaxy from the IR SED modelling. 

IC 883 being less luminous in the IR, but farther away than 
Arp 299-A and in a much more advanced merger stage, contains 
both a starburst and an AGN. The nuclear starburst in IC 883 
is app arently less prolific than in Arp 299-A (iPerez-Torres et all 
2009). However, the AGN in IC883 seems to be stronger than 
the one in Arp 299-A (IPerez-Torres et a l. 2010), and even pow- 
ering the radio luminosity at circumnuclear scales, although its 
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contribution to the emission of the entire galaxy at IR wave- 
lengths is apparently insignificant. 

4. Summary 

We have imaged the nuclear and circumnuclear regions of the 
LIRG IC 883 at radio frequencies. Our e-EVN observations, to- 
gether with archival VLBI observations, reveal the presence of 
at least four non-thermal components in the nuclear regions of 
IC883. Three of these are transient sources and one is a long- 
lived, variable compact source, which is very likely an AGN. 
Our e-MERLIN observations highlight the presence of a striking 
double-sided structure, likely representing a warped disk/ring, in 
agreement with previous studies on the interstellar medium of 
this galaxy. 

The source we identify as an AGN candidate is powering the 
radio emission in IC883 at nuclear and circumnuclear scales, 
and yet, as indicated by our radio and NIR observations, as 
well as by our model of the IR SED, this LIRG displays very 
active star formation both within the innermost lOOpc of its 
nucleus, and in the circumnuclear regions. We have estimated 
lower limits for the CCSN rate in the innermost nuclear re- 
gions of 0.6^Q3yr _1 based on the detection of three radio SNe 
above a 5<x luminosity threshold of 4 x 10 27 ergs~' Hz -1 , and 
1-1-0 6 y r 1 f° r ^ e ent i re g a l ax y. based on the nuclear radio SNe 
and the NIR SNe 2010cu and 201 1 hi discovered in the circum- 
nuclear regions. Our most recent NIR Gemini-North observa- 
tions of IC 883 clarify the nature of SN 201 lhi as a Type IIP SN, 
in agreement with its non-detection in our observations at radio 
frequencies. 

The composite starburst-AGN nature of IC 883, revea led 
previously by optical spectroscopy (e.g.. lYuan et al.lfeOlOl) . is 
now supported by our radio observations for the first time, thus 
making IC 883 one of the few galaxies where the starburst-AGN 
connection can be studied in detail. 
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